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Full Diversity Unitary Precoded Integer-Forcing 

Amin Sakzad and Emanuele Viterbo^^ 


Abstract —We consider a point-to-point flat-fading MEMO 
channel with channel state information known hoth at transmit¬ 
ter and receiver. At the transmitter side, a lattice coding scheme 
is employed at each antenna to map information symbols to 
independent lattice codewords drawn from the same codebook. 
Each lattice codeword is then multiplied by a unitary precoding 
matrix P and sent through the channel. At the receiver side, 
an integer-forcing (IF) linear receiver is employed. We denote 
this scheme as unitary precoded integer-forcing (UPIF). We show 
that UPIF can achieve full-diversity under a constraint based 
on the shortest vector of a lattice generated by the precoding 
matrix P. This constraint and a simpler version of that provide 
design criteria for two types of full-diversity UPIF. Type I uses 
a unitary precoder that adapts at each channel realization. Type 
II uses a unitary precoder, which remains fixed for all channel 
realizations. We then verify our results by computer simulations 
in 2 X 2, and 4x4 MIMO using different QAM constellations. 
We finally show that the proposed Type II UPIF outperform the 
MIMO precoding X-codes at high data rates. 

Index Terms —MIMO, Integer-Forcing, unitary precoding, lat¬ 
tice codes, full-diversity. 

I. Introduction 

Multiple-input multiple-output (MIMO) channels with mul¬ 
tiple antennas appeared in the early 1990’s as the key technol¬ 
ogy to achieve high spectral efficiencies in wireless channels. 
These channels are known to obtain higher capacity gains 
in comparison to single antenna point-to-point channels. For 
a MIMO channel matrix with n* transmit antennas 

and rir receive antennas, a diversity gain of at most ntrir 
is achievable. If an encoder/decoder pair gives the maximum 
possible diversity then the entire system is said to 

achieve full-diversity. Let rig denote the number of symbols 
transmitted per channel use. If Ug = mm{nr,nt}, then the 
encoding scheme is said to have full-rate. 

For MIMO channels with channel state information at the 
receiver (CSIR), carefully designed space-time block-codes 
(STBC) can provide full-diversity m. Algebraic full-diversity 
full-rate STBCs with large coding gain were later designed 
in la, a. The design criteria for these codes are known 
as the rank and the non-vanishing determinant criteria IIJ, 
0. In order to take full advantage from these codes, an 
optimal maximum likelihood (ML) detector such as a sphere 
decoder a, with high computational complexity should be 
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implemented at the receiver side. Other well-known linear 
receivers 0, including zero-forcing (ZF) and minimum mean 
square estimation (MMSE) linear receivers, trade-off error 
performance for reduced computational complexity. Specific 
families of full-diversity STBCs were investigated and con¬ 
structed in 0-ill] for use with linear receivers (ZF and 
MMSE). In these schemes, full-diversity is obtained at the 
cost of limited rate. 

Precoding schemes that assume CSIR, channel state infor¬ 
mation at the transmitter (CSIT), and an ML decoder usually 
make use of a singular value decomposition (SVD) to convert 
the MIMO channel into parallel subchannels. Then appropriate 
precoders employed over these parallel channels to minimize 
the error probability of the system over different constellation 
sizes. Some examples are the well-known full-diversity E- 
dmin precoders for 4-QAM 112, and MIMO precoding X- 
and Y- codes Ga, which are known to provide a diversity 
order of = {nr — ^ -f l) (n* — ^ -f l), and the 

full-diversity precoders introduced and investigated in m. 
Note that dnt,nr.,ns shows a trade-off between full-rate and 
full-diversity. In one hand, if rit = Ur and we use a full- 
rate MIMO precoding X- or Y- code, then the diversity order 
is dnj,n,.,ris = + l) < n^ (full-diversity). On the other 

hand, if we want to ensure full-diversity for MIMO precoding 
X- or Y- codes, we need Ug = 2, which means the system 
has very limited rate. Eor a detailed comparison between the 
ML-decoding complexity of the above three schemes, we refer 
the reader to Table III of llT4l . 

On the other hand, in the presence of CIST and having a 
linear receiver such as MMSE linear receiver at the destina¬ 
tion, optimal precoder designs have been extensively studied 
( ifTSll - ifTSl ). The design criteria for these schemes include 
maximizing the rate or minimizing the error performance. 
The MIMO linear precoding techniques such as regularized 
ZE mi are alternative approaches that can provide full- 
diversity in conjunction with linear receivers, see ll20l and 
references therein. The decoding complexity of such schemes 
is relatively low in comparison to ML-decoding algorithms but 
incur some error performance degradation. 

In this paper, assuming CSIT and CSIR we propose a 
full-diversity unitary precoding scheme for n x n MIMO 
(ur = rit = n), where the receiver is equipped with the 
recently proposed integer-forcing (IE) linear receiver El- 
In IE framework, the transmitter employs a lattice coding 
scheme, and sends independent lattice codewords from the 
same codebook, across different transmit antennas (or layers). 
At each receive antenna, we recover an integer linear com¬ 
bination of the lattice codewords depending on the entries 
of the channel matrix H. As all the lattice codewords are 
drawn from the same lattice codebook, the integer linear 
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combinations are themselves guaranteed to be lattice points. 
To further estimate the original transmitted lattice codewords, 
a system of linear equations with integer coefficient matrix A 
should be solved at the receiver. In fact, the IF linear receiver 
exploits the linearity of lattices to hrst eliminate the noise 
and then bareness the interference between transmitted data 
across different layers. The key step underpinning IF is the 
selection of A to approximate the channel matrix H. Hence, 
a matrix B is needed such that A is full rank and BH « A, 
with minimum quantization error at high signal-to-noise ratio 
(SNR) values. The problem of hnding A based on H using 
lattice reduction algorithms is addressed in ll22]l . 12^ . The 
integer-forcing linear receiver is also shown to provide a full 
receive diversity order Ur = n and a full multiplexing gain 

in im, EH. 

Using perfect STBCs like the Golden code ||25]| . ||3l . the 
authors of 1261 have shown that under an IF decoder, a constant 
gap to capacity is attainable. A perfect STBC code is a lin¬ 
ear dispersion space-time code over a Quadrature-Amplitude 
Modulation (QAM) constellation, which is full-rate, satisfies 
non-vanishing determinant property, and its generator matrix 
is a unitary matrix. In this paper, we aim at relaxing the hrst 
two constraints of perfect STBCs and keeping the last one 
only. We analyze the diversity order of a lattice space-time 
encoding scheme, which is precoded by a unitary matrix and 
uses an IF linear receiver. First, singular value decomposition 
(SVD) precoding is used to diagonalize the channel, and then 
an optimal precoder matrix is designed. Our optimization 
criterion is based on maximizing the minimum distance of 
the input lattice, similar to the approach of 1271 . We note that 
differently from l27l . our precoder matrices are unitary and 
our decoder is the IF linear receiver and not the ML receiver. 
Related works focus on a different design criterion based on 
maximizing the mutual information ( l28l . l29l . l30l . and OTl 
and references therein). 

From the receiver point of view, using IF linear receiver 
instead of ML decoders in slow-fading channels is advanta¬ 
geous in terms of complexity. For the IF receiver, the lattice 
reduction algorithm is performed only at the beginning of 
each quasi-static interval for which the channel H is assumed 
to remain constant over multiple codeword transmissions. 
Subsequently, for each codeword transmission, only a system 
of linear equations has to be solved to recover an estimate 
of that codeword. On the other hand, for the ML decoding, 
a sphere decoder should be used for each codeword within a 
quasi-static channel interval. 

We summarize the contributions of this paper as follows; 

• Unitary precoded integer-forcing (UPIF) - Considering 
the availability of CSIT, the SVD decomposition WSV^ 
for H„xn is first performed. Using V at the transmitter 
and W at the receiver, the channel is reduced to a diag¬ 
onal channel matrix S. An additional unitary precoding 
matrix P, is then employed as VP. At the destination, 
the receiver is equipped with an IF linear receiver. 

• Full-diversity constraint - We find an upper bound 
on the probability of error of UPIF (Theorem [T]) and 
show analytically that if P satisfies a specific constraint 
(Theorem |^, the full diversity (n^) can be achieved by 


our UPIF. The exact constraint is hard to verify, since it 
is related to finding the minimum distance of a lattice 
depending on S and P. Theorem provides a simpler 
constraint, based on the non-vanishing minimum product 
distance of the lattice generated by P and guarantees the 
full-diversity of UPIF. 

• Precoder design criteria - The constraints on P provide 
the design criteria for full-diversity-achieving precoders 
for IF MIMO linear receivers. We define two types of 
full-diversity UPIF schemes: (i) the precoder matrix P 
depending on the channel matrix singular values (diago¬ 
nal elements of S) is adapted at each quasi-static channel 
interval {Type I UPIF), (ii) a fixed unitary precoder P 
(independent of S) is used for all channel realizations 
(Type II UPIF). For Type I UPIF, finding an analytical 
optimization technique for P remains an open problem. 
For Type II UPIF we propose the use of full-diversity 
lattice generator matrices with non-vanishing minimum 
product distance El-illTl. 

• Simulation results - Monte-Carlo simulations are pre¬ 
sented in support the derived theoretical results. For 
example, it is shown that a 2 x 2 full-diversity algebraic 
rotation matrix given in llJTl in a Type II UPIF scheme 
with a 64-QAM constellation, can perform as close as 
IdB away from ML decoder for the same code and 2dB 
better than MIMO precoding X-codes ifTSll . 

Notation. The superscripts (•)^ and (•)^ denote transposi¬ 
tion and Hermitian transposition, respectively. Let Q and Q' be 
a group and its subgroup respectively, then Q/Q' denotes the 
quotient group. The sets Z, C, K, and Z[i] denote the ring of 
rational integers, the field of complex numbers, the field of real 
numbers, and the ring of Gaussian integers, respectively, where 
= —1. Real and imaginary parts of a complex number 
z are denoted by 3fl(z) and 3(z), and \z\ and arg(z) are 
the modulus and the unique phase, respectively. The notation 
||v|| stands for the Euclidean norm of a vector v £ C”. A 
k X k matrix X = [x^,... ,x^]^ is formed by stacking the 
fc—dimensional row vectors Xi,...,Xfc, and denotes the 
k X k identity matrix. Finally, let v be a vector, then its j-th 
entry is represented by [v]j. 


II. System Model 

We first review the notion of complex lattices lf3^ which 
are essential for the rest of the paper. 

Definition 1: Given a set of d-dimensional vectors 
{li,... ,ik] in C'^, a fc-dimensional complex lattice A is the 
set of points 

^ ^ X'm^m ■ X.ra £ ^ £ C /■ . 

m—1 J 

More briefly we can write A = {zL: z £ Z[z]^}, where 
L = [£j,..., represents a basis of the lattice A. In 

other words, every point x £ A can be represented as a 
Gaussian integer linear combination of basis vectors. For any 
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point X S A the Voronoi cell V(x) is 
V = ^ amt-m - llv-x|l < ||v-y|j, Vy S A, e c| 

m—1 } 

The Voronoi cell for the origin is denoted by V. Let L be a 
matrix with as its rows, then L is called the generator 
matrix of the lattice Al- Throughout the paper, we only 
consider full rank lattices where d = k. For example 
is a lattice with the identity matrix as generator matrix. Let 
S C C", then the linear span of S is defined as: 

span(5) = < ^ CpVp: q gN, Vp G S, Cp G C 

[p=i 

Let Br{0) = {x G ||x|| < r}, be the d-dimensional ball 
of radius r centered at 0 and dim(V) denotes the dimension 
of a subspace V of For a d-dimensional lattice Al, we 
define the m-th successive minima, for 1 < to < d as 

em(AL) = inf {r: dim (span (Al n 0 ^( 0 ))) > to} . ( 1 ) 


The TO-th successive minima of Al is the infimum radius r 
such that there are to independent vectors of Al in Sr( 0 ). 

A d-dimensional lattice Al is called full-diversity if for 
all distinct x,y G Al, [x]m ^ [y\m for all 1 < to < d. 
The minimum product distance of a full-diversity lattice Al is 
denoted by dp,min(AL) (l35l): 

dp,min(AL) = min TT|[x]m|. (2) 

xGAlMO} 

A lattice code C C A is a finite set of points of A. One 
way of constructing lattice codes is through sublattices (nested 
lattices). A subset A' C A is called a sublattice if A' is a lattice 
itself. Given a sublattice A', we define the quotient A/A' as 
a lattice code. The fine lattice A is called the coding lattice 
and the coarse lattice A' is usually referred as the shaping 
lattice. The lattice code C = A/A' corresponds to a finite 
constellation of lattice points carved from the lattice A, within 
a fixed Voronoi cell of A'. In other words, let V' be the Voronoi 
cell of origin for lattice A', then C = A n V'. Therefore, the 
shape of such constellation is governed by V'. For example, 
the Voronoi region of the 4-dimensional checkerboard lattice 

A' = D 4 = I (a:i, • • • ,Xi) gZ^: ^ Xi = even 

t 771—1 

has exactly 24 faces, with the relevant vectors being the 
24 minimum Euclidian lattice points 1391. This lattice is a 
sublattice of 4-dimensional integer lattice, i.e. D 4 C and 
hence so is the scaled lattice 8D4 = { 8 x: x G D4} c z^. 
The lattice code A/A' = I 8 D 4 includes all 4-dimensional 
integer points inside the Voronoi region of the lattice 8D4. 

In general, a common choice for the sublattice A' of A is 
the scaled lattice gA for some integer p G Z as in 1381 . As 
another simpler example, for A = Z[i]‘^ and A' = gZ[i]‘^, 
we have A/A' = Zg[i] for which Zg is the ring of integers 
modulo g. For g = 2 the cubic shaped lattice code Z 2 [i] is 
simply a 4-QAM constellation. For the rest of this paper we 
use A = Z[i]‘^ and choose g to be a power of 2. 


We consider a quasi-static flat-fading n x n MIMO chan¬ 
nel, where the channel state information is available at both 
transmitter and receiver. The channel matrix is denoted by 
H G C"^", where the entries of H are i.i.d. complex Gaussian 
random variables ^ CAf{0, 1). We use an n-layer lattice cod¬ 
ing scheme, where the information transmitted across different 
antennas are independent. For 1 < m < n, the m-th layer 
is equipped with a lattice encoder E : Ef C" which 
maps a message G over the ring TZ into a lattice 
codeword x™ G A/A' C C" in the complex space. The 
matrix X = [x}’,..., x^]^ is actually a space-time codeword, 
where its rows are all lattice codewords. This lattice space- 
time codeword will then be precoded using a unitary matrix. 
The precoding matrix U can be derived using the components 
of the channel matrix H and another unitary matrix P, to 
be optimized later. Together with U, the matrix X forms a 
space-time codeword UX to be sent through the channel. 

Let H = WSV^ be the singular value decomposition 
(SVD) of the channel matrix where W,V G are 

two unitary matrices and S is a diagonal matrix given by 
S = diag(cfi,..., (fn) with cfi > • • • > (t„ all in M. A unitary 
precoder matrix 

U = VP (3) 


is then employed at the transmitter where P G is a uni¬ 

tary matrix that needs to be optimized. If X = [xf,..., x^]^ 
denotes the matrix of transmitted vectors, the received signal 
Y is given by 

Y = Vp ■ HUX + Z, (4) 

SNR 

where p = denotes the average signal-to-noise ratio at 
each receive antenna and the entries of Z are i.i.d. distributed 
as CA/^(0,1). Let M be the generator matrix of a complex 
lattice Aj^, then by using the standard conversion from 
complex lattice to the equivalent real one, we have the real 
lattice generator matrix as: 


/ 5i(M) S(M) \ 

-3(M) SR(M) ) ■ 

With the above transformation, 0 can be written as 


Y = • HUX -f Z, (5) 

where all the matrices Y,H, U,X, and Z are in K 2 nx 2 n ^ 
suitable block diagram is then as Fig. [T] Upon receiving Y 
at the destination, we multiply it by to get Y' = W^Y. 
Substituting U from Q into Q the channel can be modeled 
as: 

Y' = Vp ■ SPX -f Z', ( 6 ) 

where Z' = W^Z. Note that Z' continues to be an i.i.d. 
complex Gaussian matrix with entries ^ CAA(0,1) because 
W is unitary. 

At this point, the receiver employs integer-forcing. That is 
a linear receiver architecture that creates an effective integer¬ 
valued channel matrix. The entries of this effective channel 
matrix are used to recover integer combinations of the code¬ 
words. This is in contrast to the other linear receivers and 
ML-based decoders, which attempt to recover the transmitted 
codewords directly. A feature of integer-forcing linear receiver 
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Fig. 1. Real Block diagram of unitary precoded integer-forcing. 


is in using the same lattice codes at both the transmitter 
as well as the receiver. In this framework, the sender uses 
a layered transmission scheme, and transmits independent 
lattice codewords simultaneously across the layers. At the 
receiver side, each layer is allowed to decode an integer 
linear combination of transmitted lattice codewords. Since 
any integer linear combination of lattice points is another 
lattice point, the decoded point will be another lattice point. 
Assuming that the effective channel matrix is full rank, these 
integer combinations can then be solved for the original lattice 
codewords. 

The goal of integer-forcing linear receiver is to project SP 
(by left multiplying it with a receiver filtering matrix B) onto 
a non-singular integer matrix A. In order to uniquely recover 
the information symbols, the matrix A must be invertible over 
the ring TZ. Thus, we have 

Y" = BY' = ^ ■ BSPX + BZ'. (7) 

For the IF receiver ETIl formulation, the model is given by 

Y" = v^-AX+ 7^-(BSP-A)X +BZ' 

'-..-^ 

quantization noise term 

^ Vp-AX-I-E, (8) 

where • AX is the desired signal component, and the 
effective noise is 


(BSP-A)X-f-BZ'. (9) 

We further denote the effective noise term along the m-th 
layer by and by stacking these row vectors, we get E. The 
average energy of the effective noise along the m-th row of 
Y" is defined as 


G(a„,bm) = pljb^SP - -f |jb„^|p, (10) 


where a^ and b^ denote the m-th row of A and B, 
respectively. Note that in order to increase the m-th layer 
effective signal-to-noise ratio lISTIl 


SNR,„ 


P 

G(a^, bjT^) 


( 11 ) 


the term G(am,bm) has to be minimized, for each m, by 
appropriately selecting the matrices A and B. We refer to the 
above signal model as unitary precoded integer-forcing. Note 
that the signal strength in the numerator of is simply 


p (instead of p||am|P) due to the fact that a^x is itself a 
codeword, which is decodable up to an effective noise variance 
of p. Specifically, for lattice codes, it is often easier to think of 
each codeword as surrounded by a ball of radius ^Jnp. If the 
noise falls within this ball, then decoding is successful. The 
received signal-strength viewpoint is more appropriate for i.i.d. 
random coding methods, where a joint typicality analysis can 
be applied. 


A. Decoding Complexity Advantage of UPIF 

We highlight the decoding complexity advantage of UPIF 
in comparison to ML decoder techniques such as sphere 
decoder H in slow-fading MIMO channels. In these channels, 
the phases and the amplitudes of the channel coefficients can 
be assumed approximately constant over a long period. One 
such period is called a quasi-static channel interval or simply 
an interval. Let the channel coefficients remain unchanged for 
a time period t, which is much larger than the codeword 
length 2 n, so that many codewords can be transmitted in 
this interval. For example, let m = ^ G Z be the number 
of feasible codeword transmissions within one such quasi¬ 
static channel interval. In IF, the matrices A and B are 
computed once, based on H and estimated p in the first 
codeword transmission. In the subsequent m — 1 codeword 
transmissions within the same interval, only a system of 
linear equations, with integer coefficients A, will be solved 
to recover the original information symbols at the receiver. 
Referring to Table II of ||22l, the overall complexity of using 
IF in an interval is dominated by the complexity of finding 
the matrix A as the m systems of 2n linear equations with 
2n unknown variables can be solved using Gauss-Jordan 
elimination process with computational complexity of order 
m X (2n)^. This is in contrast to ML decoding algorithms in 
slow-fading MIMO channels. In these scenarios, an algorithm 
such as sphere decoder should be employed for each and 
all m codeword transmissions inside one quasi-static channel 
interval. Therefore, the overall computational complexity for 
an interval is m times the complexity of the sphere-decoder 
algorithm ||4|. This makes IF linear receiver computationally 
more efficient than ML decoding algorithms in slow-fading 
channels. 
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III. Diversity Analysis 

We first recall the definition of diversity and a weak version 
of Woodbury identity HOl . which are used in the rest of this 
Section. 

Definition 2: In an 2n x 2n MIMO system and at a high 
signal-to-noise ration p, if the average probability that 
a transmitted symbol vector is wrongly decoded, is approx¬ 
imated by {c.p)~^, then 5 is called the diversity order (or 
diversity gain) and c is called the coding gain. For a MIMO 
system with precoding, if (5 = ( 2 n)^, then, we say that the 
precoder achieves full-diversity order. 

Lemma 1 (Woodbury Identity): Let Mi and M 2 be two 
d X d invertible matrices, then the following identity holds: 

Id - Ml (Id + M2 Mi)-^ M 2 = (Id + MiM2)-^ . 

We now proceed to analyze the diversity order of UPIF. The 
effective noise in is not Gaussian distributed due to the 
quantization noise term. Note that matrix ( 12 ^ + p ■ 
is a positive definite matrix as it is a diagonal matrix with 
positive entries, hence it admits a Cholesky decomposition 


(l2„ + p-S'*S) ^=LL\ 

(12) 

The optimum value of that minimizes 

is m 

(ED given a„ 

bjn = P • {^2n + P ' 5jP (SP) 

(13) 

By defining 

S4l2„ + p.SP(SP)\ 


we re-write (pj]) as 


= p-a„SP'*S-i. 

(14) 

Let 

Lp ^ P'^L, 

(15) 

then the average energy of the effective noise term along the 


m-th layer is given by 


p||b„SP-a,„f+ ||b,^f 

p-a„(l2„-(SP)'‘S-iSP)a(; 

(16) 

P ■ Wn {^2n + P ■ (SP) SP^ a 

™(17) 

p-a„P'‘(l2„+p-S'*S) 


(18) 


(19) 


where (ED holds because of ( pAl i and equation (9) in fT2\ . 
(ED is true because of the Woodbury identity in Lemma E] for 

d = 2n, Mi = ^.(SP)\ M 2 = Vp-SP, 

and (El and (El are obtained based on (ED and (ED, 
respectively. We denote the probability of error for decoding 
the m-th layer in the integer lattice by (m, SP,Z^"). 
It follows that 

Pe (m,SP,Z2-) = Pr (^|e^| > , (20) 


where is the m-th row of E defined in (H. 

The following Theorem provides an upper bound on 
Pe(m,SP,Z2-). 

Theorem 1 (Upper Bound on Probability of Error): For 

all 1 < m < 2 n, the term Pe (m, SP,Z^”) as in ( | 20 | , is 
upper bounded as 

Pe (m,SP,Z 2 ") < exp (-ce^„_^+i(AL-i)) , ( 21 ) 

where c is a constant independent of p and e 2 n-m+i(^L “0 
is the {2n — m + l)-th successive minima of the lattice with 
generator matrix defined in ( Ep . 

Proof: The proof is given in Appendix Note that in 
the proof we assume the use of dithering which makes x„j 
a random variable distributed uniformly within the shaping 
region. ■ 

We are interested in Pg (2n, SP, Z^") which is an upper 
bound for Pg (m, SP,Z^"), for all 1 < m < 2n. We define 
the error probability for unitary precoded integer forcing over 
Z^" as 

Pg (SP, Z^") 4 Pg (2n, SP, Z^”) . (22) 

Based on ( | 2 T] ), we have 

Pg(SP,Z) < exp (-ce?(AL-i)) . (23) 

Let the average probability 

Pg=EH (Pg(SP,Z)), 

where the expectation is taken over all channel matrices H. 

Theorem 2: Let the precoding matrix P be such that 
[Pv]i f 0, where v S Z^" is the vector satisfying ef(AL-i) = 
||Lp ^v|p, then the unitary precoded integer-forcing explained 
through (H to ED achieves full-diversity {2n)‘^. 

Proof: The proof is given in Appendix [B] ■ 

Based on the proof of the above theorem, to achieve full- 
diversity, all we need is to select a unitary precoder matrix 
P such that if V G Z^" be the vector, which arises to the 
minimum distance of the lattice A^-i, then it satisfies [Pv]i f 
0. However, enforcing this condition seems to be unpractical 
because it involves computation of dmin which is an 

NP-hard problem. For that reason, we choose to work with a 
stronger condition that is dp_i„in(Ap) f 0. It is clear that if 
<^p,min(Ap) f 0 then [Pv]j^ f 0 for every vector 0 7 ^ v € 
Z 2 ”. 

Theorem 3: Let the precoding matrix P be such that 
dp,min(Ap) f 0 , then the achievable diversity of the unitary 
precoded integer-forcing explained through (ID to ED is 

i2n)f 

According to the above theorem, if lattice codes along with 
a non-zero minimum product distance precoder are used at 
the transmitter of a MIMO channel and integer-forcing is 
employed at the receiver, full (transmit and receive) diversity 
can be achieved. We next provide two different design criteria 
for obtaining unitary precoders suitable for integer-forcing 
linear receivers. 
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IV. Optimal Design of Full-Diversity Unitary 
Precoders 

Based on the above two theorems, we next design classes 
of full-diversity unitary precoders for integer-forcing linear 
receivers. We hrst propose a method of designing optimal 
precoders based on Theorem and then we provide our 
optimal precoders based on Theorem 


A. Design of Type I UPIF 

Based on ( |2T| and Theorem]^ the optimal Type I UPIF is 
as follows: 


Pi,„p. = arg max min ||L ^Pv|p, (24) 
P 6 C' 2 „ vGZ2"\{0} 

[Pv]i#0 


where C>2n is the orthogonal group of all 2n x 2n orthogonal 
matrices with matrix multiplication operation. In other words, 
we should design a precoder matrix P such that the minimum 
distance of the lattice A^-i is maximized, while [Pv]i f 0 
for V e and ef(AL-i) = ||L“^v|p. 

Remark 1: If we relax the constraint on P € C>2n to be 
Tr (PP"^) < pq, where Tr(.) denotes the trace of a matrix, 
and remove the constriant [Pv]i f 0, we get a more general 
optimization problem, which is the main subject of study 
in na. The optimal solution for this problem is derived to be 
the hexagonal lattice 1431 for 2 x 2 real MIMO channels 
and Schlafli lattice I?4 l43l for 2 x 2 complex MIMO channels. 
However, the generator matrix of these lattices are not unitary. 
In this paper, we focus on solving the optimization problem 
in ( |24| ), which is slightly different from that of l42ll . 

The problem stated in ( |24l l seems to be a hard combinatorial 
problem, which we could not solve analytically. However, a 
systematic approach of hnding P,„p, is to perform an exhaus¬ 
tive search in Oin- A convenient parametrization of C>2n using 
(2n)^ parameters is presented in l44l . Such a parametriza¬ 
tion uses the fact that an orthogonal transformation P can 
be composed from elementary orthogonal transformations in 
two-dimensional sub-spaces. The induced matrices using the 
algorithm given in l44ll will be uniformly distributed with 
respect to Haar measure. 

For the 2x2 MIMO system, since we imposed the orthogo¬ 
nality constraint on P, we parameterize it using a single angle 
Q as follows: 


P{0) = 


cos 9 sin 0 \ 
— sin 9 cos 9 I ’ 


(25) 


It can be easily seen by the symmetry of integers that it 
suffices to consider 9 G [0,7r/4] for the maximization in ( [24l i. 
Hence, we numerically search for 

Pl5=arg max min ||L-ip(0)vf, (26) 

p(e)G02 [P(e)v]i#o 

where O 2 is the group of all 2 x 2 orthogonal matrices with 
matrix multiplication operation. This exhaustive search can 
be done by discretizing 9 G [0, 7r/4] using fine steps of 
0.001 radians and consequently hnding the minimum distance 
of Ay-i using Gauss reduction algorithm ||45]| . The 2x2 

■L'p 

real precoding matrix, which results in the highest minimum 



Fig. 2. The variation of tan0 based on the variation of tan 77 in a 2 x 2 
complex MIMO Channel using real Type I UPIF. 


distance of the corresponding lattice A^-i is considered as a 

2x2 real approximately optimal percoder matrix P[*] to be 
employed for both real and complex 2 x 2 UPIF. Using ([T^, 
the matrix is given by 




L-^ = 


paf 


0 


v/i" 






cos r] 
0 


0 

sin?7 


where ^ = 


^/2~+~^c 


af) and 


^ (v/l + P^i'j + (v/l + P'^2 


T] = tan 


-1 


yi + 




,yr+ 


per? 


It follows that L ^ is equal to 


/ cos 7 ] 

0 ) 

f cos 9 

sinP \ 

1 0 

sin 77 J 

1 — sin 9 

cos 9 J 


L-^P(P)=e 


In Fig. 1^ we plot the variation of tan 9 based on tan rj for 
10® samples of 2 x 2 complex MIMO systems employing real 
2x2 precoders P(0). In particular. 


tanp = 








paf 




where /3 l-i is the condition number of For large values 
of p, tanp ~ ^ where /3 h is the condition number of 

the channel H. The generator matrix of the lattice A^-i can 
be considered now as the equivalent channel matrix. It can be 
seen that if the condition number is less than >/3, then the best 
precoder obtains by putting 9 = 7r/4. In Fig. we plot the 
variation of the coding gain of the lattice A^-i versus tanp. 
The coding gain formula is: 


tUl-i) = 




(27) 


det (L^ 
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Fig. 3. The variation of 7 (A, - 1 ) based on the variation of tan 77 in a 2 x 2 
complex MIMO Channel using real Type I UPIF. 



Fig. 4. The histogram of 7 (Aj^-i) in a 2 x 2 complex MIMO Channel 
using real Type I UPIF. 


The coding gain measures the increase in density of A^-i over 
the integer lattice with 7 (Z^”) = 1. Based on Fig. we 
can safely say that the equivalent precoded channel matrix 
Lp ^ is at least perform as good as an integer 2x2 channel 
matrix. In other words, the precoder matrix will change the 
bad channels to equivalent channels as good as an integer 
one. Fig. |4| shows the histogram of 7 (Al-i) for the same 
10® channm samples. There is a probability mass at around 
1, which states that unitary matrices has cleared all the bad 
channels with 7 (Ah) < 1 to convert them to an equivalent 
channel matrix with 7 (Al-i) > 1. Note also that the 
maximum attainable coding gain (Hermite’s constant El) 
for 2-dimensional lattices is As it is clear from Figs, 
and 1^ this quantity has been achieved by few lattices in our 
simulations too. 


B. Design of Type II UPIF 

In previous Subsection, we introduced Piop, as a suitable 
precoder matrix for UPIF. Equations ( |25] l and ([26] l show that 
Type I UPIF precoder matrices Pi^p, and pf^fare adapted 
based on and consequently H and p. In this Subsection, 
however we introduce Type II UPIF, where the optimal pre¬ 
coder matrix Pjop, is fixed and it does not change by varying 
the channel matrix H. 

By further expanding ei(AL-i) in ( |23| ), we get: 


e?(AL-0 = \\Lf\f 

= v%-%-1v = v"(LpU)-A 
= v'‘P'‘((l2„ + p-S'*S)”^)”'pv (28) 
= v'‘PUl2„ + p-S'*S)Pv, (29) 


where ( [28] l follows from ( [T9| ). Since l 2 „-l-p S^S is a positive 
definite matrix, we can perform a Cholesky decomposition to 
obtain l 2 „ -f p ■ = DD^. With this, ( |29| ) can be written 

as ||DPv|p. Hence, we get 


e?(AL-) = 


IIDPvf 

2n 

E [Pv]: 


> 2n n n (30) 

\m—l m—1 / 

1 

( 2n 2n \ 2 n 

n n 

m—1 m—1 / 

2n 

2n n «n.in(Ap))’^ (31) 


> 


m—1 

2np(det(H))" (Aj 


(32) 


where ( [30l l follows from Cauchy-Schwarz inequality, and © 
is true because 1 + p(T^(H) > p(T^(H), for all 1 < to < 2n 
and (Ap) < [Pu]^ for all u G Z^". 

Based on ( [32l l and since almost surely det (H) f 0, the 
optimal Type II UPIF can be obtained by solving the following 
optimization problem: 


P 2 .„p, = arg ^max d;^-^ (Ap), (33) 

The solution for the maximization in ( |33] l is provided in 1321- 
13^ using algebraic number theoretic lattices. A list of full- 
diversity algebraic rotations is available in 1371 . All the lattices 
provided there satisfy non-vanishing minimum product dis¬ 
tance criterion too. Hence, thanks to l|37l, the optimal Type II 
UPIF along with their corresponding dp min are available. We 
now proceed to verify this result by conducting simulations. 


V. Simulation Results 

The UPIF scheme and MIMO precoding X-codes and Y- 
codes ifTTll share similar properties, which make them suitable 
for comparison: (i) both schemes use SVD decomposition 
technique to transform the channel matrix into a diagonal 











Fig. 5. Type I UPIF in comparison with, X-Precoders decoded with sphere 
decoding algorithm, and Type II UPML in a 2 x 2 complex MIMO Channel. 


one, (ii) the precoder matrices in both systems must be 
unitary/orthogonal matrices, (Hi) both the detectors at the 
receiver side, i.e. lattice reduction based IF linear receiver ll22ll 
and a combination of two 2-dimensional ML decoders, provide 
full receive diversity in 2 x 2 MIMO. 

For comparison we show the performance of UPIF precod¬ 
ing when the IF receiver is replaced by an ML decoder. We 
will denote this with the acronym UPML. In the following 
subsections, we simulate and compare our proposed Type FII 
UPIF decoded using lattice reduction based IF receiver, with 
Type I/II UPML, and MIMO X-Precoders/X-Codes. 

A. Type I UPIF versus X-Precoders 

We show simulations for a 2 x 2 MIMO UPIF channel over 
4/16/64-QAM constellations. The matrices A and B for IF 
linear receiver were found using Algorithm2 presented in 1i22\ . 
We have chosen our Type I UPIF by running an exhaustive 
search based on the design criterion given in ( |2^ . For ref¬ 
erence purposes, we also presented MIMO X-precoder ifTSll 
scheme decoded under sphere decoding algorithm a. The 
step size for our brute force search for both the cases (Type 
I UPIF and X-precoders) is 0.001 radians. Fig. shows the 
codeword error rate fT2\ curves versus signal-to-noise ratio 
p of the different schemes. Since the performance of UPIF 
decreases by growing p parallel to ML decoded curve, full- 
diversity of UPIF is guaranteed. In particular, at codeword 
error rate 10“^ over 4/16/64-QAM constellations, the gaps 
of 2.2dB, 2.15dB, and 2.1dB between Type I UPIF and X- 
precoder are observed. This shows a decrease in coding loss 
by increasing the constellation size. The 2x2 X-precoders 
outperform the other schemes because of two main reasons; 
(/) these precoders are designed for the particular lattice 
constellation (for example, 4-QAM, 16-QAM, or 64-QAM) 
while Type I UPIF are designed for the (infinite) lattice 
(ii) since we are only able to simulate Type I UPIF over a 
2x2 MIMO channel, the diversity gain of Type I UPIF over 



Fig. 6. Type II UPIF in comparison with X-Codes and Type II UPML in a 
2x2 complex MIMO Channel. 


X-precoders cannot be captured. This will be more visible in 
the next subsection, where Type II UPIF outperform X-codes 
in MIMO channels with larger number of antennas employed 
at larger constellation sizes. 

B. Type II UPIF versus X-Codes 

We have conducted simulations for 2x2 and 4x4 MIMO 
UPIF channel over 4/16/64-QAM constellations. For 2x2 
MIMO channels, the matrices A and B for IF linear receiver 
were found using Algorithm2 presented in ll22ll while for larger 
MIMOs we have employed Algorithm 1. We have chosen our 
Type II UPIFs based on the design criterion given in ( [3^ . The 
2x2 and 4x4 full-diversity algebraic rotation EH matrices 
which has been used for simulations have the following 
minimum product distances 0.668740, 0.438993, and 

0.289520. For reference purposes, we also presented MIMO 
X-codes US decoded under sphere decoding algorithm 0. 
For 4-QAM constellation the optimal angle for a X-Code is 
26.6 degrees while for 16-QAM it is equal to 15 degrees 113. 
We also found the optimal angle for 64-QAM to be equal 
to 8 degrees. In Figs. and |7] we plot the codeword error 
performance of Type II UPIF versus MIMO X-codes and Type 
II unitary precoders decoded using a maximum likelihood 
decoder (Type II UPML) such as sphere decoding algorithm 
versus signal-to-noise ratio Es/Nq = p. As it is clear from 
Figs. 1^ and |7j the Type II UPIF outperform MIMO X-codes 
for larger constellation sizes in higher MIMO dimensions. The 
main reason is that the Type II UPIF precoder matrices are 
designed to achieve full-diversity and the diversity gain shows 
itself when the number of antennas increases. In particular, 
the diversity order na of the optimal MIMO X-codes (in 
terms of diversity not the error performance) is known to be 
dnt,n^,n, = ^ + l) (tt* “ + l), where Ur and n* 

denote the number of receive and transmit antennas, respec¬ 
tively, and Us denotes the number of symbols. In our setting. 
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Fig. 7. Type II UPIF in comparison with X-Codes and Type II UPML in a Fig. 8. Type I versus Type II UPIF and UPML schemes in a 2 X 2 complex 
4x4 complex MIMO Channel. MIMO Channel. 


for 2x2 and 4x4 MIMO X-codes, we have ^ 2 , 2,2 = 4 and 
<^ 4 , 4,4 = 9 in comparison with fall-diversity orders 4 and 16 
for Type II UPIF, respectively. 

C. Type I UPIF versus Type II UPIF 

In this subsection, we compare Type I and Type II UPIF 
together. Fig. shows the codeword error performance of 
Type I and Type II UPIF in comparison to their corresponding 
UPML. It is obvious that the Type II UPIF and UPML 
outperform Type I UPIF and UPML precoders by no more 
than 0.5dB and IdB, respectively. 

The sub-optimality of unitary precoder matrices obtained 
using quantized exhaustive search algorithm is the main reason 
for having Type II UPIF performing better than Type I UPIF. 
Furthermore, we note that the design criteria given in ( |26l l and 
( |3^ are different and we cannot say which one can perform 
better than the other by just looking at these equations. Hence, 
considering the simplicity of Type II UPIF in comparison to 
Type I and their better error performance, we suggest using 
of Type II UPIF for practical purposes. 

VI. Summary and Directions for Future Work 

A unitary precoding scheme (UPIF) has been introduced for 
transmission over a flat-fading MIMO channel with both CSIT 
and CSIR, and an IF linear receiver. The diversity gains of 
the proposed approach have been analyzed both theoretically 
and numerically. Two different design criteria for two types 
of unitary precoders. Type I UPIF and Type II UPIF, were 
given. For both cases, an SVD is performed first to transform 
the channel to a diagonal one. Type I UPIF includes precoder 
matrices adapting based on the channel matrix singular values, 
while Type II UPIF consists of precoding matrices, which 
remain fixed by changing the channel matrix. 

Designing full-diversity unitary precoders or more generally 
space-time block codes for IF receivers without CSIT would 
be of interest. Another direction is to let the transmitter 


have access to limited feedback over a delay-free link from 
the IF receiver ED, El. In this case, designing a suitable 
codebook of unitary precoding matrices attaining higher rates 
and obtaining higher coding gains in error performance seems 
to be a promising research topic. 


Appendix A 
Proof of Theorem[T] 

For the proof of this we follow the lines of the proof of 
Theorem 7 in El- Since the minimum Euclidean distance 
of Z is unity, an error is declared if Therefore, 

Pe (m, equals to 


=Pr 



2Pr 




(34) 


<2 min 

t>0 


=2 min 

i>0 


E(exp(tem)) 

exp(^) 

E(exp (t^'(b„ 


(35) 




exp 


(^) 


,(36) 


where follows from the symmetry of effective noise 
around zero, follows from Chernoff’s bound. Since 
and are independent, the equation ( (36] t can be written as 


min 

t>o 


E(exp(ty/p ■ (b^SP-a„t, x,„)))E(exp (t ■ (b^, z'^})) 
4 exp (4*) 


(37) 

We continue to use moment generating functions of uniform 
and Gaussian distributions to further upper bound the two 
exponential terms in ( |J7l . Since the components of are 
distributed based on a Gaussian distribution with zero mean 
and unit variance, we have 


E(exp(f • (b^,z^))) < exp ^ (38) 

for the second term in We now investigate its 

first term. Let = ty/p ■ (b^SP — a^). The 
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therefore 


entries of x™ G are uniformly distributed over 

the hypercube [-f^||q„||,f^pllq^li] 

E (exp(f^ • (q^,Xm))) equals to 

= E exp ■ y^[qm]j[x„]j 


i=i 


2n 


< 


< 


J| E (exp {ty/p ■ [q„]j [x„]j)) 
i=i 

|| sinh (fVP|[qm]i[x™]j|) 


t^p\ [q: 


< exp 


i=i 

2n 

J|exp 

tVllqmf 


mlj I 


6 


(39) 

(40) 

(41) 

(42) 


where ( [39| ) follows from the independence of components of 
Xm, ( |40| l follows from the moment generating function of a 
uniform distribution, and ( |4T] l follows from the inequality 

sinh(j/) 

y 

for every y G K. Combining ( |42| i and ( |38| l in we get 
Pe(m, SP,Z) less than or equal to 


< exp ( — 


< 2 min • 

i>0 


exp 


t^p\\hm^P-aLrr 


exp 




exp 


exp 

2 


( 41 ) 


2 mm - 

exp (^) 

(43) 

t>0 

2 exp ^ 

4G(a„,b™)) ’ 

(44) 


where ( |44l i obtained by optimizing t, that is f = - 


2G(a 


Now since G{am, hm) needs to be minimized, and 
should be chosen appropriately as in m and Ea, thus we 
get 

G{arn, bm) 2 / A A 
- - - — 

where e^(ALh) denotes the m-th successive minimum of the 
lattice 

Ala = {GPd I VdGZ^"}. 

Here L~^ is a generator of the dual lattice A* a. Thus we have 
the relation (see Lemma 4 in 11211 1 

(2n)^ + (3n)^ (2n)^ + (3n)^ 


i(Alj) < 


'=-2n-m+lV'-L 


(A^) 


'■2n—m+1 


(Al-)’ 


(45) 


where h) is the (2n—m+l)-th successive minima 

of the lattice A^^. Therefore, we have 


P 


bj^) 


> 


'-2n—771+1 


(Al-) 


Co 


( 46 ) 


where cq = (2n)^ + (3n)^. Using inequality of ( |46l l in ( |44l l, 
the probability of error for decoding the m-th layer is upper 
bounded as 

Pe {m, SP,Z2”) < exp (-ce^„_„+i(AL~i)) , (47) 

where c = +. This completes the proof 


Appendix B 
Proof of Theorem|2] 

Toward the proof of this theorem, we make use of the 
following two lemmas. The first lemma gives the expression 
for the first-order expansion of the marginal probability density 
function (PDF) of eigenvalues of H. 

Lemma 2: Let the entries of the matrix H be i.i.d. complex 
Gaussian random variables with zero mean and unit variance. 
The first-order expansion of the marginal PDF of the fc-th 
largest eigenvalue Afe of the complex central Wishart matrix 
HH^ is given by 

Fx,{\k) = qk\t +o(xt), 

as Afe 0, with dk = (2n — k + 1)^ — 1 and qk being positive 
constants. 

Lemma 3: For a scalar channel modeled by y = VSNR^x-f 
n , where n ^ A/'(0,1), and E[|xp] = 1 and a = /3'^ 
is a nonnegative random variable whose probability density 
function (PDF) Fa{a) is such that as 

Fa{a) = qa* + o(a*), (48) 

the average word error probability (WEP) Pe, which is given 
by 

pOO 

Pe = E {Pg{a)) = / exp (—fcaSNR) Fa{a)da, 

Jo 

is such that as SNR —> 00 

Pe = -f o (^SNR"(*+A^ _ 

Now we proceed to give the proof of Theorem Let v be 
the Gaussian vector for which 

= l|L; + f 

= v+'‘((l2„ + p-S''S)”^)”'pv (49) 
= (l2„ + p • Pv 

4 IIDPvf > [B]l, [Pv]J (50) 

= (1 + p+(H)) [Pv]? > pa?(H) [Pv]J(51) 
> y+MH) [Pv]^ (52) 

where follows from follows from Cholesky 

decomposition of l 2 „ + p ■ and the fact that D is a 

diagonal matrix, ( |5T] i is true because 1 -f pa^ (H) > p^im 
and ( [ 5 ^ follows from the fact that crJ(H) = Ai(H). 


^?(Al-0 
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Since the quantity [Pv]i is not zero. Using Lemmas 

and ( |5^ , we have 

Pe = Eh (Pe(SP, E^")) < + O , 

where to is a function of n, qi dehned in ( |48l ), and [Pv] i 7 ^ 0. 
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